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We present the theory of inclusive decays of hadrons containing a heavy quark and discuss its most important 
applications to the decays of B mesons. We also review the theoretical understanding of the hadronic parameters 
Ai and A2 (or fi^ and /ig) entering the heavy-quark expansion. 



1. INTRODUCTION 

Hadronic bound states of a heavy quark with 
light constituents (quarks, antiquarks and gluons) 
are characterized by a large separation of mass 
scales: the heavy-quark mass mg is much larger 
than the mass scale Aqcd associated with the 
light degrees of freedom. Equivalently, the Comp- 
ton wave length of the heavy quark (Aq ~ I/tuq) 
is much smaller than the size of the hadron con- 
taining it (i?had ~ 1/Aqcd)- In such a situa- 
tion, it is appropriate to separate the physics as- 
sociated with these two scales, so that all depen- 
dence on the heavy-quark mass becomes explicit. 
The framework in which to perform this separa- 
tion is the heavy-quark (or l/mq) expansion [y- 
|8) , which is a specific realization of the operator 
product expansion (OPE) ||, |o). 

There are (at least) two important reasons why 
it is desirable to separate short- and long-distance 
phenomena: A technical reason is that after this 
separation we can calculate a big portion of the 
relevant physics (i.e. all short-distance effects) 
using perturbation theory and renormalization- 
group techniques; in particular, in this way we are 
able to control all logarithmic dependence on the 
heavy-quark mass. An important physical reason 
is that, after the short-distance physics has been 
separated, the long-distance physics may simplify 
due to the realization of approximate symmetries, 
which relate the long-distance properties of many 
observables to a small number of hadronic matrix 
elements. The second point is particularly excit- 
ing, since it allows us to make statements beyond 
the range of applicability of perturbation theory. 



In our case, an approximate spin-flavour sym- 
metry is realized in systems in which a single 
heavy quark interacts with light degrees of free- 
dom by the exchange of soft gluons ||ll[-|l^ . The 
origin of this symmetry is not difficult to under- 
stand. In a heavy hadron, the heavy quark is 
surrounded by a most complicated, strongly inter- 
acting cloud of light quarks, antiquarks and glu- 
ons. However, the fact that the Compton wave- 
length of the heavy quark is much smaller than 
the size of the hadron leads to simplifications. 
To resolve the quantum numbers of the heavy 
quark would require a hard probe; the soft glu- 
ons exchanged between the heavy quark and the 
light constituents can only resolve distances much 
larger than l/niQ. Therefore, the light degrees of 
freedom are blind to the flavour (mass) and spin 
orientation of the heavy quark. They experience 
only its colour field, which extends over large dis- 
tances because of confinement. In the rest frame 
of the heavy quark, it is in fact only the electric 
colour field that is important; relativistic effects 
such as colour magnetism vanish as mq —> 00. 
Since the heavy-quark spin participates in inter- 
actions only through such relativistic effects, it 
decouples. It follows that, in the limit mg — > 00, 
hadronic systems which differ only in the flavour 
or spin quantum numbers of their heavy quark 
have the same configuration of their light degrees 
of freedom. Although this observation still does 
not allow us to calculate what this configuration 
is, it provides relations between the properties of 
such particles as the heavy mesons B, D, B* and 
D* , or the heavy baryons Ab and Ac. 



Heavy-quark symmetry is an approximate sym- 
metry, and corrections arise since the quark 
masses are not infinitely heavy. Nevertheless, re- 
sults derived on the basis of heavy-quark symme- 
try are model-independent consequences of QCD 
in a well-defined limit. The symmetry- breaking 
corrections can, at least in principle, be studied 
in a systematic way. A convenient framework for 
analysing these corrections is given by the heavy- 



quark effective theory (HQET) |16|, which pro- 
vides a systematic expansion around the limit 
niQ — > GO. In the HQET, a heavy quark in- 
side a hadron moving with velocity v is described 
by a velocity-dependent field hy subject to the 
constraint '^ hy = hy. This field is related to 
the original heavy-quark field by a phase redef- 
inition, so that it carries the "residual momen- 
tum" k = pq — mgv, which characterizes the in- 
teractions of the heavy quark with gluons. The 
effective Lagrangian of the HQET is ||l6|-||l|] 
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where D'^ = 13'' — {v ■ D) u'' contains the com- 
ponents of the gauge-covariant derivative orthog- 
onal to the velocity, and gsC" = i[D'',D''] is 
the gluon field-strength tensor. The leading term 
in the effective Lagrangian, which gives rise to 
the Feynman rules of the HQET, is invariant 
under a global SU{2nh) spin-flavour symmetry 
group, where n^ is the number of heavy-quark 
flavours. This symmetry is explicitly broken by 
the higher-dimensional operators arising at order 
1/toq, whose origin is most transparent in the 
rest frame of the heavy hadron: the first oper- 
ator corresponds to (minus) the kinetic energy 
resulting from the motion of the heavy quark in- 
side the hadron (in the rest frame, («-D_l)^ is the 
operator for — k^), and the second operator de- 
scribes the chromo-magnetic interaction of the 
heavy-quark spin with the gluon field. The coef- 
ficients Ckin and Cinag result from short-distance 
effects and, in general, depend on the scale at 
which the operators are renormalized. At the tree 

level, Ckin = Cmag = 1- 

At this point it is instructive to recall a more 



familiar example of how approximate symmetries 
relate the long-distance properties of several ob- 
servables. The strong interactions of pions are 
severely constrained by the approximate chiral 
symmetry of QCD. In a certain kinematic regime, 
where the momenta of the pions are much less 
than IGeV (the scale of chiral-symmetry break- 
ing), the long-distance physics of scattering am- 
plitudes is encoded in a few "reduced matrix el- 
ements", such as the pion decay constant. An 
effective low-energy theory called chiral pertur- 
bation theory provides a systematic expansion of 
scattering amplitudes in powers of the pion mo- 
menta, and thus helps to derive the relations be- 
tween different scattering amplitudes imposed by 
chiral symmetry |20|] . A similar situation holds 
for the case of heavy quarks. Heavy-quark sym- 
metry implies that, in the limit where niQ ^ 
Aqcd , the long-distance properties of several ob- 
servables is encoded in few hadronic parameters, 
which can be defined in terms of operator matrix 
elements in the HQET. 

In particular, the forward matrix elements of 
the two dimension- five operators in (hi). 



Okin = -hy{iDA_fhy 
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play a most significant role in many applications 
of the HQET. They appear, for instance, in the 
spectroscopy of heavy hadrons and in the descrip- 
tion of inclusive weak decays M . For the B me- 
son, we define two parameters Ai and A2 (or, 
equivalently, /ij and /i^) by |2l| 



- Ai = /i^ = 



1 
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Whereas A2 is directly related to the mass split- 
ting between vector and pseudoscalar mesons 
through 



m%, — 171% ~ 4A2 + 0(l/mb) 



(4) 



the parameter Ai cannot be determined from 
hadron spectroscopy. It is, however, related to 
the difference of the pole masses of two heavy 



quarks through 

mb-rric = (jub - tjid) 
1 



Ai 



2mt 



2mD 



(5) 



Here ms = |("t-b + Sms.) and mo — |(f«_D + 
3to£i*) denote the spin-averaged meson masses, 
defined such that they do not receive a contri- 
bution from the chromo-magnetic interaction. In 
the foUowing section, we discuss the theoretical 
status of the parameters Ai and A2 and their 
properties under renormahzation. 

2. OPERATOR RENORMALIZATION 
AND THE VIRIAL THEOREM 

In quantum field theory, local composite oper- 
ators such as Okin and Omag require renormal- 
ization to be well defined. The matrix elements 
of these operators (i.e. Ai and A2) depend on 
the subtraction scheme and on the renormahza- 
tion scale fi introduced in the process of remov- 
ing ultraviolet divergences. This dependence is 
cancelled by an opposite scale dependence of the 
Wilson coefficient functions Ckin and Cmag in the 
effective Lagrangian (|l|) . The most common regu- 
larization scheme in QCD is dimensional regular- 
ization [2^, ^ with modified minimal subtrac- 
tion (MS scheme ||2J]). Ultraviolet divergences 
are regulated by working in d = 4 — 2e space-time 
dimensions, and are subtracted by removing the 
poles in 1/e. The scale dependence of the Wilson 
coefficients is governed by renormalization-group 
equations of the form 



^^^c{^l)^^ci^l)., 



(6) 



where 7 is called the anomalous dimension. 

The anomalous dimension of the kinetic oper- 
ator Okin vanishes to all orders in perturbation 
theory, 7kin = 0. This is a consequence of the 
reparametrization invariance of the HQET, i.e. 
and invariance under infinitesimal changes of the 
velocity v used in the construction of the effective 
Lagrangian pa-[E7|. It follows that in dimen- 
sional regularization Ckin = 1- The anomalous 
dimension of the chromo-magnetic operator does 
not vanish, however. It has been known at the 



one- loop order since many years |T^, Q , but only 
very recently the two-loop coefficient has been 
calculated |E8, 29 . For N„ = 3 colours, the result 



7mag 
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where Uf is the number of light-quark flavours. 
Given this result, together with the one- loop 
matching condition [|7| 
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one can work out the next-to-leading order ex- 
pression for Cmag ifJ-) 1 starting from the general 
solution to M) m,m: 



sM 



C(^) = C{mQ) exp 



da. 



aairriQ) 



7(Q^s) 
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where l3{as) = dasin) / din /j, is the /3 function. 
The product of the Wilson coefficient C,nag(A*) 
and the scale-dependent matrix element A2 (p-) is 
renormalization-group invariant. Hence, in the 
presence of renormahzation effects, what can be 
determined from the vector-pseudoscalar mass 
splitting in M) is the combination 



A2 = C„,ag(mb) A2(mfc) « 0.12 GeV' 



(10) 



The numerical value is specific for the B system, 
since according to our definition the renormalized 
parameter A2 depends logarithmically on the b- 
quark mass. 

Unfortunately, the question about the value of 
the kinetic-energy parameter Ai is more difficult 
to answer, even from a conceptual point of view. 
We have already mentioned that spectroscopic re- 
lations involving Ai, such as (^, depend on the 
heavy-quark pole masses, which are not physi- 
cal parameters. In the remainder of this section, 
we shall argue that indeed the parameter Ai is 
not physical, but has an intrinsic ambiguity of 
order Aq(-,j3. (However, the difference of the val- 
ues of Ai in two different hadrons is a physical pa- 
rameter and can be extracted from spectroscopy.) 
This means that Ai must be defined in a non- 
perturbative way; in other words, different defini- 
tions cannot be simply related to each other using 



Table 1 

Theoretical estimates of the parameter Ai 
(QCDSR: QCD sum rules, HQSR: heavy-quark 
sum rules, Exp.: experimental data on inclusive 
decays, QM: quark models) 



Reference 


Method 


-Al [GeV2] 


Eletsky, Shuryak Q 
BaU, Braun |||l 
Neubert || 


QCDSR 
QCDSR 
QCDSR 


0.18 ±0.06 
0.52 ±0.12 
0.10 ±0.05 


Gimenez et al. [|4 


Lattice 


-0.09 ±0.14 


Bigi et al. m 


HQSR 


> 0.36 


Gremm et al. [|6) 
Falk et al. ^ 
Chernyak [|38[ 


Exp. 
Exp. 
Exp. 


0.19 ±0.10 
«0.1 
0.14 ±0.03 


Hwang et al. Q 
Be Fazio [|o| 


QM 
QM 


0.5 ±0.1 
0.66 ±0.13 



perturbation theory. The wide spread in the the- 
oretical predictions for the parameter Ai, shown 
in Tab. |^, is partially a reflection of this prob- 
lem. Future efforts should concentrate on under- 
standing better the relations between the various 
definitions underlying these estimates. 

The problem of the ambiguity in the value of 
the heavy-quark kinetic energy is closely related 
to the well-known ambiguity in the definition of 
the pole mass of a heavy quark |41, 4^]. The rea- 
son lies in the divergent behaviour of perturbative 
expansions in large orders, which is associated 
with the existence of singularities along the real 
axis in the Borel plane, the so-called renormalons 
Q-[^. When the pole mass rriQ is related to 
a short-distance mass (such as the running mass 
in the MS scheme), which is a well-defined quan- 
tity in perturbation theory, then the correspond- 
ing perturbation series 

"^Q = "^Q° {l + cias{mQ) + C2al(mQ) + .. . 

+ c„<(mQ) + ...}, (11) 

contains numerical coefficients c„ that grow as 
n\ for large n, rendering the series divergent and 



not Borel summable. The best one can achieve is 
to truncate the series at the minimal term, but 
this leads to an unavoidable arbitrariness of or- 
der Atuq '^ Aqcd (the size of the minimal term) . 
This observation, which at first sight seems a se- 
rious problem, should not come as a surprise. We 
know that because of confinement quarks do not 
appear as physical states in nature. Hence, there 
is no unique way to define their on-shell proper- 
tics such as a pole mass. In view of this, it is 
actually remarkable that QCD perturbation the- 
ory "knows" about its incompleteness and indi- 
cates, through the appearance of renormalon sin- 
gularities, the presence of non-perturbative ef- 
fects. We must first specify a scheme how to 
truncate the QCD perturbation series before non- 
perturbative parameters such as mq and Ai be- 
come well-defined quantities. The actual values 
of these parameters will depend on this scheme. 

In the difference of the pole masses on the left- 
hand side in (|^), the leading renormalon ambi- 
guity of order Aqcd cancels. However, there re- 
mains a residual ambiguity of order Aqcq/toq 
Ip^ . This point was unclear for some time, since 
the corresponding renormalon singularity does 
not appear in the so-called bubble approxima- 
tion, in which most explicit calculations in the 
Borel plane are performed [^ [5^. However, it is 
now known that this singularity must be present 
beyond the bubble approximation, and the ambi- 
guity in the difference mi, — rric is then absorbed 
by a corresponding ambiguity of order Aq^j-, in 
the parameter Ai in (||) Q. 

It is instructive to consider the same problem 
from a different point of view, by studying the 
properties of the kinetic operator under renormal- 
ization in regularization schemes with a dimen- 
sionful cutoff parameter /x. In such schemes, the 
operator Okin can mix with the lower-dimensional 
operator hyh^ (the "identity"), because the two 
operators have the same quantum numbers. If 
such a mixing is present, it leads to an ad- 
ditive contribution to the parameter Ai of the 
form fp C[as{fJ,)], which one would like to sub- 
tract in order to define a renormalized parame- 
ter. The coefficient C can be calculated order 
by order in an expansion in the small coupling 
constant as(/i), and it appears at first sight that 



the quadratically divergent term could be sub- 
tracted using perturbation theory. This impres- 
sion is erroneous, however, because C may con- 
tain non-perturbative contributions of the form 
exp[-87r//3oas(^)] = (Aqcd/m)^. which cannot 
be controUed in perturbation theory [ p4[ . Such 
terms contribute an amount of order Aqqq to the 
parameter Ai, which is of the same order as the 
renormahzed parameter itself. Hence, if the ki- 
netic operator mixes with the identity, it is nec- 
essary that the quadratically divergent contribu- 
tion to Ai be subtracted in a non-perturbative 
way, and hence the heavy-quark kinetic energy 
by itself is not directly a physical quantity. 

The question whether there is a mixing of the 
kinetic energy with the identity, and whether 
there exists a corresponding renormalon singu- 
larity, has been addressed by several authors, 
with seemingly controversial conclusions. At the 
one-loop order, such a mixing has indeed been 
observed when the HQET is regularized on a 
space-time lattice iQ. Likewise, a "physical" 
definition of a parameter Ai(/i) has been sug- 
gested, which absorbs certain 0{as) corrections 
appearing in the zero-recoil sum rules for heavy- 
quark transitions pq] . This definition is such that 
dAi(^)/d/^^ ex as(/z), indicating again a one- loop 
mixing of the kinetic energy with the identity. 
On the other hand, this mixing has not been 
observed at the one-loop order in two Lorentz- 
invariant cutoff regularization schemes |p3] , which 
use a Pauli-Villars regulator or a cutoff on the 
virtuality of the gluon in one-loop Feynman dia- 
grams 1^^. This observation appeared as a puz- 
zle, which was resolved when it was shown that 
the mixing between the kinetic energy and the 
identity is forbidden at the one-loop order in all 
regularization schemes with a Lorentz-invariant 
UV cutoff,[| but that in general there is no sym- 
metry that protects the matrix elements of the 
kinetic operator from quadratic divergences, and 
so a mixing with the identity occurs from the two- 
loop order on B2]. 



^The fact that a mixing at the one-loop order was observed 
in Refs. | p5| , M is a consequence of the exphcit breaking of 
Lorentz invariance by the regularization schemes adopted 
in these calculations. 



In understanding this situation, the virial the- 
orem of the HQET is of great help. This theorem 
provides a relation between the kinetic energy of 
a heavy quark inside a hadron and its chromo- 
electric interactions with gluons [^. For the B 
meson, it can be written in the form 



lim 



l\2 _ 1 



3 



[v ■ v') 

{B\ Okin \B) 



(12) 



Note that in the rest frame of the initial or 
the final meson, the operator appearing on the 
left-hand side of this relations only contains the 
chromo-electric field El. = — G"*. The virial 
theorem is a most useful relation. Not only is 
it intuitive, generalizing a well-known concept 
of classical physics; it also helps in understand- 
ing better the properties of the kinetic operator. 
For instance, it has been used to estimate the 
hadronic parameter Ai using QCD sum rules [p3|. 
More important, however, is the fact that the 
virial theorem can be employed to analyse the 
properties of the kinetic operator under renor- 
malization p2|. (Parenthetically, we note that 
this theorem has also been used in the calcu- 
lation of the two-loop anomalous dimension of 
the chromo-magnetic operator [ESl M.) In the 
limit v' —> V, the properties of the two op- 
erators Okin a-nd O^^ = hyiigsG'^'^hy are re- 
lated to each other. The chromo-electric operator 
O'^^ can mix with the lower-dimensional operator 
iv^v'^ —v'^v'^) hy'hy. However, in any regulariza- 
tion scheme with a Lorentz-invariant ultraviolet 
cutoff, such a mixing can only come from Feyn- 
man diagrams involving gluons attached to both 
heavy-quark lines, since otherwise there is no way 
to get the factor (v'^v"^ — v'^v''^). Such diagrams 
appear first at the two-loop order, as shown in 
Fig. |l|. The virial theorem thus provides for a 
simple explanation of the fact that the mixing of 
the kinetic operator with the identity was not ob- 
served at the one-loop order in Lorentz-invariant 
regularization schemes. On the other hand, an 
explicit calculation confirms that the mixing is 
present at the two-loop order, and thus there is 
indeed a quadratic divergence (and thus a renor- 
malon problem) affecting the matrix elements of 




Figure 1. Two-loop diagrams contributing to 
the mixing of the chromo-electric operator with 
lower-dimensional operators. Not shown are two 
copies of the first two diagrams with the "outer" 
gluon attached to the other heavy-quark line. 



the kinetic operator p2[ . 

At the end of this discussion, we stress that 
the "renormalon ambiguities" are not a concep- 
tual problem for the heavy-quark expansion. It 
can be shown quite generally that these ambigu- 
ities cancel in all predictions for physical quanti- 
ties [^-|Q. The way the cancellations occur is 
intricate, however. The generic structure of the 
heavy-quark expansion for an observable is of the 
form: 



Observable '^ C[as{mQ)] I 1 



A 



, (13) 



where C[as{mQ)] represents a perturbative co- 
efficient function, and A is a dimensionful non- 
perturbative parameter. The truncation of the 
perturbation series defining the coefficient func- 
tion leads to an arbitrariness of order Aqcd/'tiq, 
which cancels against a corresponding arbitrari- 
ness of order Aqcd in the definition of the non- 
perturbative parameter A. Thus, only when the 
short- and long-distance contributions are com- 
bined in the heavy-quark expansion, an unam- 
biguous result is obtained. 

3. INCLUSIVE B DECAY RATES 

Inclusive decay rates determine the probability 
of the decay of a particle into the sum of all pos- 
sible final states with a given set of global quan- 
tum numbers. An example is provided by the 
inclusive semileptonic decay rate of the B meson, 
T{B ^ X^ii'), where the final state consists of 
a lepton-neutrino pair accompanied by any num- 
ber of hadrons with total charm-quark number 
rir = 1. Here we shall discuss the theoretical 



description of inclusive decays of B mesons (an 
analogous description holds for all hadrons con- 
taining a heavy quark) |6l|-||7C)|l. From the the- 



oretical point of view, such decays have two ad- 
vantages: first, bound-state effects related to the 
initial state (such as the "Fermi motion" of the 
heavy quark inside the hadron Bq, p9J) can be ac- 
counted for in a systematic way using the heavy- 
quark expansion; secondly, the fact that the fi- 
nal state consists of a sum over many hadronic 
channels eliminates bound-state effects related to 
the properties of individual hadrons. This sec- 
ond feature is based on a hypothesis known as 
quark-hadron duality, which is an important con- 
cept in QCD phenomenology. The assumption 
of duality is that cross sections and decay rates, 
which are defined in the physical region (i.e. the 
region of time-like momenta), are calculable in 
QCD after a "smearing" or "averaging" proce- 
dure has been applied [ [7l| . In semileptonic de- 
cays, it is the integration over the lepton and 
neutrino phase space that provides a "smearing" 
over the invariant hadronic mass of the final state 
(so-called global duality). For non-leptonic de- 
cays, on the other hand, the total hadronic mass 
is fixed, and it is only the fact that one sums over 
many hadronic states that provides an "averag- 
ing" (so-called local duality). Clearly, local dual- 
ity is a stronger assumption than global duality. 
It is important to stress that quark-hadron du- 
ality cannot yet be derived from first principles, 
although it is a necessary assumption for many 
applications of QCD. The validity of global du- 
ality has been tested experimentally using data 
on hadronic r decays 172]. Some more formal at- 
tempts to address the problem of quark-hadron 
duality can be found in Refs. |7^, |74| . 

Using the optical theorem, the inclusive decay 
width of B mesons can be written in the form 

r(B ^X) = —ImlBlTlB), (14) 

where the transition operator T is given by 



d''xT{C,s{x),C,s{0)} 



(15) 



In fact, inserting a complete set of states inside 
the time-ordered product, we recover the stan- 
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bgsT^u^Gf'b 



Figure 2. Perturbative contributions to the tran- 
sition operator T (left), and the corresponding 
operators in the OPE (right). The open squares 
represent a four-ferniion interaction of the effec- 
tive Lagrangian £cff, while the black circles rep- 
resent local operators in the OPE. 
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for the decay rate. Here CcS is the effective 
weak Lagrangian corrected for short-distance ef- 
fects arising from the exchange of gluons with vir- 
tualities between raw a-nd irib |75|] - [|79| . If some 
quantum numbers of the final states X are spec- 
ified, the sum over intermediate states is to be 
restricted appropriately. In the case of the in- 
clusive semileptonic decay rate, for instance, the 
sum would include only those states X containing 
a lepton-neutrino pair. 

In perturbation theory, some contributions to 
the transition operator are given by the two-loop 
diagrams shown on the left-hand side in Fig. |[ 
Because of the large mass of the h quark, the 
momenta flowing through the internal propaga- 
tor lines are large. It is thus possible to construct 
an OPE for the transition operator, in which T is 
represented as a series of local operators contain- 
ing the heavy-quark fields. The operator with 
the lowest dimension, d = 3, is hb. It arises by 
contracting the internal lines of the first diagram. 



The only gauge-invariant operator with dimen- 
sion 4 is bilp b; however, the equations of motion 
imply that between physical states this operator 
can be replaced by mi,bb. The first operator that 
is different from bb has dimension 5 and contains 
the gluon field. It is given by bgsCr^„G^''b. This 
operator arises from diagrams in which a gluon 
is emitted from one of the internal lines, such 
as the second diagram shown in Fig. 0. For di- 
mensional reasons, the matrix elements of such 
higher-dimensional operators are suppressed by 
inverse powers of the heavy-quark mass. 

In the next step, the hadronic matrix elements 
of the local operators in the OPE are expanded 
in powers of 1/mb, using the technology of the 
HQET. The result is M, M, || 



{B\bb\B) 
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6A2 + Oil/mb) . (17) 



Thus, any inclusive decay rate can be written in 
the form il-fHil 



T{B ^ Xf) 



Gln4 
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where the prefactor arises naturally from the 
loop integrations, and c{^ are calculable coeffi- 
cient functions (which also contain the relevant 
CKM matrix elements) depending on the quan- 
tum numbers / of the final state. It is instruc- 
tive to understand the appearance of the kinetic- 
energy contribution Ai /2ml ■ This is nothing but 
the field-theory analogue of the Lorentz factor 
(1 - v^)^/^ ~ 1 — k;^/2m^, in accordance with 
the fact that the lifetime, r = l/F, for a mov- 
ing particle (the b quark) increases due to time 
dilation. 

The main result of the heavy-quark expansion 
for inclusive decay rates is the observation that 
the free quark decay (i.e. the parton model) pro- 
vides the first term in a systematic l/rrib expan- 
sion |pl|. For dimensional reasons, the corre- 
sponding rate is proportional to the fifth power of 



the 5-quark mass. The non-perturbative correc- 
tions, which arise from bound-state effects inside 
the B meson, are suppressed by two powers of the 
heavy-quark mass, i.e. they are of relative order 
(AqcD/^JT-b)^- Note that the absence of first-order 
power corrections is a consequence of the equa- 
tions of motion, as there is no independent gauge- 
invariant operator of dimension 4 that could ap- 
pear in the OPE. The fact that bound-state ef- 
fects in inclusive decays are strongly suppressed 
explains a posteriori the success of the parton 
model in describing such processes pO, pl| . 

The hadronic parameters Ai and X2 appearing 
in the heavy-quark expansion (nq) have been dis- 
cussed in the previous section. For a given inclu- 
sive decay channel, what remains to be calculated 
is the coefficient functions c^. This can be done 
using perturbation theory. We shall now discuss 
three important applications of this general for- 
malism. 

3.1. Determination of \Vcb\ from Inclusive 
Semileptonic Decays 

The extraction of |Vcb| from the inclusive 
semileptonic decay rate of B mesons is based 
on the general expressio n (Iq ), with the short- 
distance coefficients 1621 -[B4 1 



„SL 



„SL 



Wr 



cb\ 



l-Sx-" 



ix^ - x^ 



-12a:* In; 



+ OK)J , 



(19) 



Here x = mc/mb, and mt and tUc are the pole 
masses of the b and c quarks, defined to a given or- 
der in perturbation theory [p2|. The 0{as) terms 
in C3 are known exactly ||83| , while only par- 
tial calculations of higher-order corrections exist 
gj, ^. The main sources of theoretical uncer- 
tainties are the dependence on the heavy-quark 
masses, unknown higher-order perturbative cor- 
rections, and the assumption of global quark- 
hadron duality. A conservative estimate of the 
total theoretical error on the extracted value of 
\Vcb\ is 5\Vcb\/\Vcb\ ~ 10% ^- Taking the result 
of Ball et al. [p5[ for the central value, and us- 
ing tb = (1.60±0.03) ps for the average i?-meson 



lifetime |B7[, we find 



\Vcb\ = (0.040 ±0.004) 



^si 



1/2 



^10.8%, 
= (40±lexp±4th)xl0-^ (20) 

In the last step, we have used BgL = (10.8±0.5)% 
for the semileptonic branching ratio of B mesons 
(see below). The value of \Vcb\ extracted from the 
inclusive semileptonic width is in excellent agree- 
ment with that obtained from the analysis of the 
exclusive decay B — » D*iD using heavy-quark 
symmetry [p8|-[|93|. This agreement is gratifying 
given the differences of the methods used, and it 
provides an indirect test of global quark-hadron 
duality. Combining the two measurements gives 
l^cfcl = 0.039 ± 0.002 le). After Kd and K^, 
this is now the third-best known entry in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. 

3.2. Semileptonic Branching Ratio for De- 
cays into T Leptons 

Semileptonic decays of B mesons into r leptons 
are of particular importance, since they are sensi- 
tive probes of physics beyond the Standard Model 
P4| . From the theoretical point of view, the ra- 
tio of the semileptonic rates (or branching ratios) 
into T leptons and electrons can be calculated re- 
liably. This ratio is independent of the factor ttt,^, 
the hadronic parameter Ai, and CKM matrix el- 
ements. To order l/ml, one finds [^^ ^ Q 

BjB^XTi?^) _ ^_^ A^ . X 

B{B^XeI>,,) -nXc,Xr)+ ^2^9[X,,Xr) 

= 0.22 ±0.02, (21) 

where / and g are calculable coefficient functions 
depending on the mass ratios Xc — rric/mh and 
Xr = nir/rnb, as well as on as{mf)). Two new 
measurements of the semileptonic branching ratio 
of b quarks, B{b -^ XtDt), have been reported 
by the ALEPH and OPAL Collaborations at LEP 
1^, ||. The weighted average is (2.68 ± 0.28)%. 
Normalizing this result to the LEP average value 
Bib^ Xeue) = (10.95 ±0.32)% |1, we obtain 



B{b^ XrDr) 



0.245 ±0.027, 



(22) 



B{b^ XeDr) 

in good agreement with the theoretical prediction 

(pl|) for B mesons. 



3.3. Semileptonic Branching Ratio and 
Charm Counting 

The semileptonic branching ratio of B mesons 
is defined as 



BsL 



V{B -^ Xeu) 



Y.^V{B^Xlv) + T 



had 



(23) 



where Fhad and Fiarc are the inclusive rates for 
hadronic and rare decays, respectively. The main 
difficulty in calculating i3sL is not in the semilep- 
tonic width, but in the non-leptonic one. As men- 
tioned previously, the calculation of non-leptonic 
decay rates in the heavy-quark expansion relies 
on the strong assumption of local quark-hadron 
duality. 

Measurements of the semileptonic branching 
ratio have been performed by various experi- 
mental groups, using both model-dependent and 
model-independent analyses. The status of the 
results is controversial, as there is a discrepancy 
between low-energy measurements performed at 
the T(4s) resonance and high-energy measure- 
ments performed at the Z^ resonance. The situa- 
tion has been reviewed recently by Richman [p^ , 
whose numbers we shall use in this section. The 
average value at low energies is i?sL = (10.23 ± 
0.39)%, whereas high-energy measurements give 
J5sl(6) = (10.95±0.32)%. The label {b) indicates 
that this value refers not to the B meson, but 
to a mixture of h hadrons (approximately 40% 
B-, 40% B°, 12% Bs, and 8% Afc). Assuming 
that the corresponding semileptonic width Fsl [b) 
is close to that of B mesons Jj we can correct 
for this fact and find Bsl — {tb/u) Bshib) = 
(11.23 ± 0.34)%, where n = (1.56 ± 0.03) ps is 
the average lifetime corresponding to the above 
mixture of b hadrons. The discrepancy between 
the low- and high-energy measurements of the 
semileptonic branching ratio is therefore larger 
than three standard deviations. If we take the 
average and inflate the error to account for this 
fact, we obtain 



BsL = (10.80 ±0.51)% 



(24) 



An important aspect in interpreting this result 
is charm counting, i.e. the measurement of the 

^Theoretically, this is expected to be a very good 
approximation. 



average number n^ of charm hadrons produced 
per B decay. Theoretically, this quantity is given 
by 

nc^l + B{B^ Xees') - B{B -^ X„oc) , (25) 

where B{B ^ Xccs') is the branching ratio for 
decays into final states containing two charm 
quarks, and B{B -^ Xnoc) ~ 0.02 is the Stan- 
dard Model branching ratio for charmless decays 
Q-llOll. The average value obtained at low en- 
ergies is jic — 1.12±0.05 IM, whereas high-energy 
measurements give n^ — 1.23 ± 0.07 |102|. The 
weighted average is 



rir = 1.16 ±0.04. 



(26) 



The naive parton model predicts that Bsl ~ 
15% and ric ~ 1.2; however, it has been known 
for some time that perturbative corrections could 
change these predictions significantly |Q. With 
the establishment of the heavy-quark expansion, 
the non-perturbative corrections to the parton 
model could be computed, and their effect turned 
out to be very small. This led Bigi et al. to con- 
clude that values Ssl < 12.5% cannot be accom- 
modated by theory |103{ . Later, Bagan et al. 
have completed the calculation of the 0(0;^) cor- 
rections including the effects of the charm-quark 
mass, findin g tha t they lower the value of Bsl 
significantly [104|. Their original analysis has re- 
cently been corrected in an erratum. Here we 
shall present the results of an independent nu- 
merical analysis using the same theoretical input 
(for a detailed discussion, see Ref. |105| ). The 
semileptonic branching ratio and ric depend on 
the quark-mass ratio rric/mh and on the ratio 
^Jb/rrib, where ii is the scale used to renormal- 
ize the coupling constant as{iJ.) and the Wilson 
coefficients appearing in the non-leptonic decay 
rate. The freedom in choosing the scale fj, reflects 
our ignorance of higher-order corrections, which 
are neglected when the perturbative expansion is 
truncated at order as. We allow the pole masses 
of the heavy quarks to vary in the range 



TO6 = (4.8 ±0.2) GeV, 
mb~ma^ (3.40 ± 0.06) GeV , 



(27) 



corresponding to 0.25 < mc/mi, < 0.33. The 
value of the difference rrif, — rric is obtained 
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from d) using Ai = -(0.4 ± 0.2) GeV^ Non- 
perturbative effects appearing at order 1/m^ in 
the heavy-quark expansion are described by the 
single parameter A2 defined in (UQ); the depen- 
dence on the parameter Ai is the same for all in- 
clusive decay rates and cancels out in the predic- 
tions for BsL and Uc- For the two choices /i — mb 
and fj, = mb/2, we obtain |105| 



B: 



SL 



12.0 ±1.0%; fi^rrib, 

10.9 ±1.0%; fi^mb/2, 

1.20 T 0.06; fi^nib, 

1.21^0.06; n^mb/2. 



(28) 



The uncertainties in the two quantities, which re- 
sult from the variation of rric/rnb in the range 
given above, are anticorrelated. Notice that the 
semileptonic branching ratio has a stronger scale 
dependence than ric- By choosing a low renor- 
malization scale, values Bsl < 12% can easily 
be accommodated. This is indeed not unnatural. 



Using the BLM scale-setting method [106], it has 
been estimated that /z > 0.32?nf, is an appropriate 
scale to use in this case [ p4| . 

The combined theoretical predictions for the 
semileptonic branching ratio and charm counting 
are shown in Fig. y. They are compared with the 
experimental results obtained from low- and high- 
energy measurements. It has been argued that 
the combination of a low semileptonic branching 
ratio and a low value of Uc would constitute a 
potential problem for the Standard Model |lOl| . 
However, with the new experimental and the- 
oretical numbers, only for the low-energy mea- 
surements a discrepancy remains between theory 
and experiment. Note that, with (25), our re- 
sults for ric can be used to calculate the branch- 
ing ratio B{B -^ Xccs'), which is accessible to a 
direct experimental determination. Our predic- 
tion of (22 ± 6)% for this branching ratio agrees 
well with the preliminary result reported by the 
CLEO Collaboration, which is B{B — > Xccs') — 
(23.9 ±3.i 



|107| 



4. SUMMARY 

We have presented the theory of inclusive de- 
cays of hadrons containing a heavy quark, and 
discussed some of its most important applications 



1.4 



1.3 



fi"l.2 



1.1 



1 
8 9 10 11 12 13 14 

BsL (%) 

Figure 3. Theoretical prediction for the semilep- 
tonic branching ratio and charm counting as 
a function of the quark-mass ratio mc/mb and 
the renormalization scale /i. The data points 
show the average experimental values obtained in 
low-energy (LE) and high-energy (HE) measure- 
ments, as discussed in the text. 
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to the decays of B mesons: the determination of 
\Vcb\ from inclusive semileptonic decays, semilep- 
tonic decays into r leptons, and the semileptonic 
branching ratio. The theoretical tools that allow 
us to perform quantitative calculations are the 
heavy-quark symmetry, the heavy-quark effective 
theory, and the I/ttiq expansion. In the case of 
inclusive decay rates, the non-perturbative infor- 
mation entering the theoretical description is en- 
coded in two hadronic parameters, Ai and A2 (or 
/ij and /ig). We have reviewed the theoretical un- 
derstanding of these parameters, concerning both 
their numerical values and their properties under 
renormalization. The parameter A2, renormalized 
at the scale /x = m(,, can be extracted from the 
spectroscopy of B mesons. The parameter Ai, on 
the other hand, suffers from a renormalon ambi- 
guity problem and thus needs a non-perturbative 
subtraction to be well defined. 
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